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ABSTRACT 

We present ground-based measurements of the transmission and emission spectra of the hot- Jupiter 
WASP-19b in nine spectroscopic channels from 1.25 to 2.35 /*m. The measurements are based on the 
combined analysis of time-series spectroscopy obtained during two complete transits and two complete 
secondary eclipses of the planet. The observations were performed with the MMIRS instrument on the 
Magellan II telescope using the technique of multi-object spectroscopy with wide slits. We compare 
the transmission and emission data to theoretical models to constrain the composition and thermal 
structure of the planet's atmosphere. Our measured transmission spectrum exhibits a scatter that 
corresponds to 1.3 scale heights of the planet's atmosphere, which is consistent with the size of spectral 
features predicted by theoretical models for a clear atmosphere. We detected the secondary eclipses 
of the planet at significances ranging from 2.2 to lAAcr. The secondary eclipse depths, and the 
significances of the detections increase towards longer wavelengths. Our measured emission spectrum 
is consistent with a 2250 K effectively isothermal 1-D model for the planet's dayside atmosphere. 
This model also matches previously published photometric measurements from the Spitzer Space 
Telescope and ground-based telescopes. These results demonstrate the important role that ground- 
based observations using multi-object spectroscopy can play in constraining the properties of exoplanet 
atmospheres, and they also emphasize the need for high-precision measurements based on observations 
of multiple transits and eclipses. 

Subject headings: planets and satellites: atmospheres — planets and satellites: individual: WASP-19b 
— techniques: photometric 



1. INTRODUCTION 

Ground-based observations using transit techniques 
have been playing an increasingly important role in the 
study of exoplanetary atmospheres over the last few 
years. It wasn't that long ago that th e first ground- 
based detections of the t ransmission (Rcdficl d et al.1 
120081: ISnellen et al.1 I2008D and thermal emission 
(Ide Mooii fc Snellen! [20091 : ISing fc Lopez-MoralesI 
l2009t iGillon et al.1 12009D of exoplanets were obtained. 
Yet today, ground-based transit measurements us- 
ing photometric techniques are now standard. The 
photometric technique has seen its widest appli- 
cation to m easurements of th e rmal emiss i on fro m 
hot-Jupiters (iRoeers et al.l [200l iCroll et al.l l2010allbl: 



Anderson et al.l 120101: iGibson et alj 120101 lAlonso et al l 



2010t iLopez-Morales et al.1 l2010t ICroll et all l2011bl: 
de Mooii et al.l 120111: ISmith et al.1 120111: iBurton et al.1 
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20ll IZhao et all I2012bl: ICrossfield et all l2012al: 
Zhao et all 12012a!: iDeming et al.l l2012t iLendl et al.l 



2013h . but it has also been used for measurements 
of transmission s pectrum of a coo l super-Earth 
(ICroll et all 120114 IdeMooii et al.1 l2012t iMurgas et all 



2012: INarita et aljl2012ir 

While photometry is useful, high-precision spec- 
troscopy over wide bandpasses is ultimately needed to 
significantly improve our understanding of exoplanet at- 
mospheres. Only these kind of data can unambiguously 
reveal atmospheric properties by resolving spectral fea- 
tures and distinguishing the overlapping lines from dif- 
ferent chemical species. Interpreting sparse photometry 
requires heavy reliance on theoretical models that have to 
make assumptions about many of the fundamental prop- 
erties of the planets that we ultimately want to determine 
observationally. Furthermore, spectroscopy opens up the 
possibility of discovering unanticipated phenomena much 
more so than photometry. 

The development of ground-based transit spectroscopy 
of exoplanet atmospheres has not enjoyed parallel matu- 
ration with the photometric technique, and such mea- 
surements are still very rare. There are two reasons 
for this. One reason is that ground-based observa- 
tions using photometric techniques can acquire simul- 
taneous observations of reference stars to enable high- 
precision relative corrections for variations in Earth's at- 
mospheric transparency, while this is not possible with 
single-object spectrographs. The second reason is that 
spectrographs typically have input feeds, either slits or 
fibers, that do not encircle all the light from a source. 
This introduces random variations in the light measured 
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TABLE 1 
Observing Log 



UT Date Event Exposures Airmass Seeing Conditions 

Time (s) # 



2012 Mar 11 00:18 
2012 Mar 13 00:01 
2012 Apr 04 00:42 
2012 Apr 05 23:23 



2012 Mar 11 06:44 Eclipse 40 

2012 Mar 13 03:55 Transit 40, 30 

2012 Apr 04 06:47 Transit 40 

2012 Apr 06 05:49 Eclipse 40 



422 1.29 -> 1.04 -+ 1.37 0.7" 

304 1.32 -> 1.04 -> 1.06 0.5" 

408 1.07 -> 1.04 -> 2.00 0.8" 

430 1.17 -> 1.04 -> 1.60 0.7" 



clear 
clear 

occasional thin cirrus 
mostly clear, passing 
cirrus for 15 min 



at the detector due to variations in seeing and imper- 
fect guiding. While there has been some development 
of post-processing algorithms to correct single-object, 
narrow-slit transit spectroscopy for variations in atmo- 
spheric transparency and slit losses (iSwain et ail 20101 : 
ICrossfield et^l201lL [20121 1 Waldmann et al.ll2012lT the 
technique has not been widely adopted. 

We have recently developed a technique for ground- 
based transit spectroscopy of exoplanet atmospheres that 
overcomes the limitations of single-object, narrow-slit ob- 
servations. The idea is to use multi-object slit spectro- 
graphs to obtain simultaneous observations of a target 
and reference stars so that differential spectroscopy anal- 
ogous to differential photometry can be performed. Light 
losses in the image plane are eliminated by using wide 
slits (> 10" is recommended). We have previously used 
this new technique to place constraints on the transmis- 
sion spectrum of t he super-Earth G J 1214 b in the optical 
and near- infrared ([Bean et al.ll20lcl l20TlT) . iGibson et all 
(|2013f) have also recently used this technique to study the 
optical transmission spectrum of the hot- Jupiter WASP- 
29. 

We present in this paper transit and secondary eclipse 
observations of the ultra short-perio d (P = 0.8d) hot- 
Jupiter WASP-19b (|Hebb et al.l(2010h . As a low-density 



hot- Jupiter (R p 



1.31 R 



1.15 M Jup ) or- 



biting very close to a star slightly cooler and smaller 
than the Sun (i?* = 0.94 i? Q ), WASP-19b is expected 
to have some of the largest spectral features in transmis- 
sion and emission among the known transiting exoplan- 
ets. WASP-19 itself is also is an intermediate-brightness 
(H = 10.6) transiting planet host star. This means that 
there are numerous similar brightness reference stars in 
the roughly 5' field-of-view typical of multi-object slit 
spectrographs. 

Observations of the transmission and emission spec- 
trum of WASP-19b have the potential to address nu- 
merous questions that have arisen concerning the nature 
of hot- Jupiter atmospheres. Our motivation was to ob- 
tain data that could be used to determine the compo- 
sition and temperature structure of the planet's atmo- 
sphere. Such information would, for example, be useful 
for constraini ng theories about the origin of thermal in- 
versions (e.g., iHubenv et all 120031: IBurrows et al.l [20071 : 
Fortnev et all 120081: iKnutson et al.l 120101 ) the statistics 
of albedo and heat redistribution (e.g ., IBurrows et al.l 



2008; IBudail 12011 ICowan fc Agdll20ll . and the possi 
bility of high carbon-to-oxygen a bundance ratios (C/ O) 
in hot- Jupiter atmospheres (e.g., Madhusudhan 2012). 

The paper is laid out as follows. We present the new 
data we have obtained for WASP-19b in §2. We de- 
scribe the analysis of the transit and eclipse light curves 



to determine the transmission and emission spectrum of 
the planet in §3. The implications of these data for the 
properties of WASP-19b's atmosphere are considered in 
§4. We conclude in §5 with a discussion of the results. 

2. DATA 
2.1. Observations 

We obtained time-series spectroscopy during two tran- 
sits and two secondary eclipses of the planet WASP- 19b 
using the MMIRS instrument ([McLeod et al.l 12012ft on 
the Magellan II (Clay) telescope at Las Campanas Ob- 
servatory in March and April 2012. A log of the obser- 
vations is given in Table [TJ We used the multi-object 
mode of MMIRS with a slit mask that allowed us to 
gather spectra simultaneously for WASP-19 and three 
other stars . As for our preyious o bservations using this 
technique ()Bean et~aL|[2010l |20iTI) , the slit widths were 
12" to avoid slit losses. The slit lengths were 30", which 
provided significant coverage of the sky background. We 
used an HK grism as the dispersive element and an HK 
filter to isolate the first order spectra. Spectra from 1.22 
to 2.36 /im with a dispersion of 6.6 A pixel -1 were ob- 
tained for all the objects. 

Complete transits and eclipses were observed without 
interruption, and a minimum of 30 minutes of data were 
obtained in each case both before ingress and after egress. 
All data were obtained at airmasses less than 2.0, with 
most obtained at airmasses less than 1.5. The expo- 
sure times used were 30 and 40 s, and the overhead was 
typically 14 s including the time to read and reset the 
detector and write the data to disk. The off-axis guide 
camera and wavefront sensor associated with the instru- 
ment were used to maintain stable pointing and the cor- 
rect shape of the primary mirror continuously during the 
observations. 

Although the wavelength range we observed in con- 
tains many narrow sky emission lines, we did not nod 
the telescope during the observations to enable pairwise 
image subtraction for removing the background. This 
staring approach is justified because the edges of the slits 
in the mask are likely reasonably parallel, and thus the 
sky lines have the same instrumental profile in the spa- 
tial direction. Furthermore, the sky emission dominates 
the background over the telescope emission at these near- 
infrared wavelengths, and thus the spatial variability can 
not be removed any better by nodding than simply in- 
terpolating over the spatially resolved background spec- 
tra obtained in each exposure. The staring approach is 
advantageous for high-precision time-series spectroscopy 
not only because it saves overhead associated with peri- 
odically moving the telescope, but also because it keeps 
the data on the same pixels throughout, and therefore 
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minimizes the influence of imperfect fiat fielding. 

The observations were obtained in generally clear con- 
ditions. There were occasional thin cirrus clouds for two 
of the nights, and drops in the fluxes can be seen in 
the raw spectrophotometry corresponding to the pass- 
ing clouds. However, these drops in fluxes are similar 
for all the stars, and the WASP-19 light curves are well- 
corrected for this effect using the reference star data. 

2.2. Data reduction 

We reduced the data using a si milar approach as f or our 
previous MMIRS observations ( Bean et al.l 1201 ID with 
two deviations. The steps in the process included col- 
lapsing the sample-up-the-ramp data cubes in to a single 
frame of total counts for the individual exposures, ap- 
plying flat field corrections based on spectroscopic flats 
taken with an internal calibration lamp, subtraction of 
the background on a wavelength-by- wavelength basis for 
each spectrum, extraction of the spectra, and determi- 
nation of the wavelength solution based on an arc lamp 
exposure taken using a mask having 0.5" wide slits. The 
deviations from our previously used data reduction algo- 
rithm are (1) ignoring the first non-destructive read in 
the up-the-ramp samples, and (2) testing the application 
of non-linearity corrections. We discuss these two issues 
in the following. 

2.2.1. Fitting the up-the-ramp samples 

The data were obtained in the "up-the-ramp" sam- 
pling mode with 5 s per read. In this mode, which is 
only possible for hybrid CMOS arrays, the detector is 
read out non-destructively multiple times during an ex- 
posure. The MMIRS system saves these images to disk, 
and so each exposure involves a sequence of full frame 
images. This sampling technique yields measurements of 
the flux as it accumulates in the pixels before they are re- 
set to begin another exposure. In the limit of a constant 
illumination and a linear pixel response, the flux values 
for a pixel will exhibit a constant rise with time over the 
course of an exposure. The advantage of this technique 
is that it can reduce the per-exposure read noise, and it 
enables checks for detector systematics. 

The first step in our data reduction is to collapse the 
corresponding data cubes of non-destructive reads to a 
single frame with the total counts per pixel for the expo- 
sure. This is done by fitting a linear trend to the ramp 
samples as a function of elapsed time on a pixel-by-pixel 
basis, and then calculating the total counts recorded for 
the pixels by multiplying the slope of this trend by the 
exposure time. 

During the analysis of the data, we noticed that the 
first read, which corresponds to a zero second integra- 
tion, always has a significantly lower bias level than the 
subsequent reads in a given exposure. Figure [1] shows 
the difference between the first two 5 s non-destructive 
reads in a dark image that illustrates this effect. Dif- 
ferences of up to and beyond 350 DN are observed, and 
the magnitude of the difference smoothly increases to- 
wards the ends of the readout channels. The dark current 
for the MMIRS detector system is estimated to be ap- 
proximately 0.01 DN s _1 pixel -1 , and is a negligible con- 
tributor to the observed difference between the first two 
ramp samples. This "bias drift" after reseting the de- 
tector is also exhibited by the HAWAII-2 RG chip in the 
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Fig. 1. — Image of the difference between the first and second 
non-destructive reads in a dark exposure. 

similar MOSFIRE instrument (jMcLean et al.ll2010f ). al- 
though the magnitude of the effect s eems to be signifi - 
cantly larger in the MMIRS detector (|Kulas et al.ll2012fl . 

The HAWAII-2 chip used in MMIRS does not have 
light-insensitive reference pixels like the later generation 
HAWAII-2 RG chips (the "R" in the RG stands for refer- 
ence pixels). Therefore, there is no straightforward way 
to calibrate the MMIRS data for this effect using exter- 
nal information. By examining dark frames, we noticed 
that the bias drift is only significant between the first 
and second non-destructive reads in dark frames, and 
the bias level is constant within the expected noise level 
for subsequent reads. Under the assumption that what- 
ever process in the detector electronics is causing this 
effect settles after < 5 s, we elected to ignore the first 
non-destructive read when fitting the up-the-ramp sam- 
ples. There is enough information with the other six or 
eight non-destructive reads per 30 or 40 s exposure, re- 
spectively, to simultaneously fit for an offset correspond- 
ing to the bias level and a slope corresponding to the 
count rate. Therefore, the final total counts recorded 
for a pixel is still the slope of the fit as a function of 
read time multiplied by the exposure time, but the zero 
point of the fit is not fixed to a known bias level like it 
otherwise could be. 

In our previous analysis of MMIRS data (|Bean et al.l 
we noticed a systematic noise pattern in the up- 
the-ramp samples in one of our data sets. The effect was 
severe enough to prompt us to consider those particular 
data unreliable, and to ignore them in our analysis. We 
note that we did not observe this same effect in the cur- 
rent data sets. All the obtained data are included in our 
analysis. 

2.2.2. Detector non-linearity 

Another deviation from our previous analysis of 
MMIRS time-series spectroscopy was to consider the in- 
fluence of the non-linear response of the detector. We 
studied the non-linear behavior of the MMIRS detec- 
tor by examining a series of flat-field exposures taken 
with the instrument in imaging mode using the inter- 
nal calibration lamp and an H filter. As is typical of 
near-infrared detectors, the MMIRS chip displays signif- 
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Fig. 2. — Top Counts in DN for non-destructive reads in a flat 
field image taken to investigate the non-linear behavior of the de- 
tector. The different point styles represent examples from the dif- 
ferent quadrants of the detector. The four examples have different 
slopes (count rates) because there is a strong instrumental through- 
put gradient across the field of view. The lines are fits to the data 
using a quadratic polynomial as a function of time (Eq. 1). The 
solid lines are the fits over the ranges where this functional form 
provides a good description of the data. The upper limits of these 
ranges correspond to the maximum values that can be corrected 
to high accuracy based these low-order fits. The data points are 
black in the well-fit regime, and grey outside. The dashed lines 
show the extrapolation of the fits. Bottom Residuals from the fit 
in terms of percentage of the data values. The point styles arc the 
same as for the top panel. 

icant non-linear response well before the full-well depth 
of the pixels is reached. The full-well is approximately 
45,000 DN (for a gain of 5 e - DN -1 ), and only 61% of 
the pixels are linear within 1% up to 16,000 DN. Fur- 
thermore, each of the four quadrants in the detector ex- 
hibits different behavior. Two quadrants even display a 
faster, rather than the expected slower, than linear rise 
in counts before turning over to plateau near saturation. 
The pixel-to-pixel variations in the non-linear behavior 
within a quadrant are also significant. 

We were aware of the approximately 16,000 DN transi- 
tion to non-linearity based on our previous observations, 
and so the exposure times of the WASP-19 data were 
tuned to keep the peak counts approximately at or be- 
low this level. Nevertheless, a significant number of pixels 
exhibit > 1% non-linearity even below this limit, and a 
few pixels near the peak in flux for the ii-band region 
of the spectra reach upwards of 20,000 DN in some ex- 
posures that were taken during periods of better than 
average seeing and low airmass. This motivated us to 
derive non-linearity corrections and apply them to the 
data to test the influence of this effect on our results. 

We determined non-linearity corrections based on a 
quadratic polynomial function. Coefficients for the cor- 
rections were derived by fitting the flux as a function of 
time for the non-destructive reads in flat field frames ob- 
tained using a total exposure time of 50 s with 1 s up-the- 
ramp samples. This was done for each pixel separately. 
The fitted function was 



where F are the fluxes of the non-destructive reads, t 
are the corresponding elapsed times, and a, b, and c are 
the parameters in the fit. The formula for correcting the 
measured fluxes is then 



F' = 



-b 2 + b\ 1- 



Ac(a-F) 



b 2 



2c 



(2) 



bt + ct 2 



(1) 



where F' is the corrected flux. 

The coefficients were determined by fitting the data af- 
ter subtracting the first non-destructive read to approxi- 
mately remove the bias. This was done because the non- 
linearity is related to the process of electron creation in 
the photosensitive layer, and this isn't influenced by the 
bias level, which arises during the reset of the detector. 
In determining the coefficients, the first non-destructive 
read wasn't fit just as for the science data, and the inter- 
cept term of the quadratic polynomial was left as a free 
parameter. 

We iteratively fit the data considering a range of pos- 
sible upper limits in count levels to determine the range 
over which the corrections were accurate using our low- 
order functional form. The median value for the upper 
limits across the detector is 32,000 DN. However, a small 
fraction of pixels are not well-described by a quadratic 
polynomial much beyond their 1% linearity ranges. Ex- 
amples of the non-linearity fits to some of the data are 
shown in Figure [2j The determined polynomial coeffi- 
cients were consistent for 50 separate exposures. The 
averages of the determined coefficients among these 50 
frames were adopted as the final correction coefficients. 

We applied the non-linearity corrections to the values 
for the non-destructive reads in each exposure before fit- 
ting for the count rate. The first non-destructive read 
was subtracted before applying the corrections to be con- 
sistent with the convention for the determination of the 
coefficients. This read was then ignored in the fit of the 
linear trend to determine the non-linear-corrected count 
rate. 

The results for the final light curves (derived tran- 
sit parameters and model fit residuals) are not signif- 
icantly different (-C lc) when using or not using the 
non-linearity corrections. This is because most of the 
pixels in the WASP-19 data sets had counts below the 
1% non-linearity level. For the results presented here, 
we did not apply the non-linearity corrections in order to 
avoid unnecessary processing steps that could add noise. 
The corrections would in principle yield more accurate 
results for more heavily exposed data, and thus could be 
of use to the wider community. Tables of the coefficients 
needed to apply Equation^ and the applicable ranges of 
the corrections are available from us upon request. 

2.3. Creation of the light curves 

An example extracted spectrum for WASP-19 from a 
40 s exposure is shown in Figure [3] This spectrum has a 
signal-to- noise ratio of 200 pixel -1 at the peak of the flux 
at 1.6 fim. We created spectrophotometric light curves 
for WASP-19 and the reference stars by summing the 
extracted spectra over wavelength. We created one, four, 
and four channels of light curves from the spectra in the 
J-, H-, and if-band atmospheric windows, respectively. 
The limits of the bandpasses in each of the atmospheric 
windows are illustrated in Figure |3l 
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Wavelength {fJ-rn) 

Fig. 3. — An example spectrum extracted from the MMIRS data 
for WASP-19. 

We divided the time-series spectrophotometric light 
curves for WASP-19 by the sum of reference star data to 
correct for the variability of Earth's atmospheric trans- 
parency during the observations. The choice of which 
reference stars to use for this correction was guided by 
which combination gave the smallest residuals in the final 
light curve fits. We found that only using the reference 
star closest on the sky to WASP-19 yielded the smallest 
model fit residuals for the J- and if-band data, while the 
sum of all three observed reference stars gave the best re- 
sults for the A"-band data. This same result was obtained 
for each of the four nights of data independently. 

The light curves for WASP-19 after correction with 
the reference star data exhibit a slow, smooth trend 
with time similar to trends that are common in ground- 
based ne ar-infra re d transit photome try (e.g.. ICroll et al.l 
I2010allbl . I2011allbl ; ISada et al.ll20Tol) . We modeled this 
trend as a second order polynomial simultaneously with 
the light curve fitting (see §3). The derived trend with 
time was not the same for the different spectral channels 
for a given event. 

We searched for additional and alternative decorrela- 
tion functions with respect to airmass and pixel position, 
but we did not detect statistically significant correlations 
with these parameters. This is in contrast t o our previ- 
ous a nalysis of MMIRS transit spectroscopy (|Bean et al.l 
1 2 1 If) . The lack of a strong correlation with airmass 
could be explained by the more similar color of WASP- 
19 with the reference stars compared to the target of our 
previous study, the mid-M dwarf GJ 1214, and its cor- 
responding reference stars. The stability of the spectral 
positions was slightly better for the current data sets, but 
not significantly so. Therefore, this likely is not the rea- 
son for the lack of a correlation with spatial pixel position 
compared to the previous study. One possible explana- 
tion for the difference is that the seeing was typically 
worse by factors of two to three for the WASP-19 data 
compared to the previously analyzed data. This means 
that the WASP-19 data were spread over more pixels 
(spatial profile FWHMs of 6 to 10 pixels compared to 
3 pixels), and this could have mitigated the influence of 
pixel-to-pixel sensitivity variations not corrected by the 
flat fielding. 

Because the light curves do not require decorrelation 
against instrument-related variables, there is little gain 
in considering data obtained significantly before or after 



TABLE 2 
Transit Parameters 



Parameter 


Value 




i (°) 
6 


78.73 ± 


0.20 


0.681 ± 


0.008 


P(d) 


0.78883910 ± 


0.00000011 


T cl (BJD TSC ) 


2455512.90256 ± 


0.00007 



TABLE 3 
Event Times 



Mid- Event Time (BJD TDB ) a 


O - C (s) 


Transits 




2455999.616301 ± 7.0E-5 


1.8±6.0 b 


2456021.703740 ±8.5E-5 


-3.0±7.3 b 


Eclipses 




2455997.6440 ± 1.0E-3 


-28 ± 86 c 


2456023.6763 ± 1.0E-3 


+21 ±89 c 



a BJDtdb is the Baryccntric Julian Date 
in the Barycentric Dynamical Time standard 
Eastman et al.ll2uTc|) . 

b Residuals from the ephemeris given in Tablc[2l 
c Residuals from the prediction of when the or- 
bital phase is 0.5 assuming a circular orbit and 
the ephemeris given in Tablc[2l and including the 
light travel time effect. 

the transits and eclipses. Indeed, the longer the extent 
of the time-series, the more likely it is that the atmo- 
spheric effect leading to the slow, smooth trend with 
time will not be well-described by a low-order polyno- 
mial. Therefore, we limited our analyses to the data 
taken within three hours of mid-transit/eclipse. The 
trimmed data still encompass more than two hours out- 
of-transit/eclipse data before and after in the cases that 
the observations spanned this length of time. 

Raw light curves for the April 5/6 secondary eclipse 
data set are shown in Figure |4j The data for WASP- 
19 are shown in this figure along with the correspond- 
ing correction functions. The correction functions are a 
combination of the reference star data, which captures 
the dominant variability, and the decorrelation against 
time. An example of the effect of passing cirrus can be 
seen by the short drop in flux about 1.5 hours after mid- 
eclipse that is common between WASP-19 and the refer- 
ence stars. 

2.4. Properties of the light curves 

The residuals from the model fits to the light curves 
when each of the two transit and secondary eclipse time- 
series are binned together at a sampling of one minute 
have rms values ranging from 650 to 1608 ppm. Most 
of the combined light curves have residuals smaller than 
1000 ppm. The residuals are 2.0 to 3.1 times the ex- 
pected level given the uncertainties estimated during the 
reduction process (including the noise contribution from 
the reference star division). In our previous analysis of 
MMIRS data, we noted that the -ff-band data had signif- 
icantly larger resdiuals than expected relative to the J- 
and A-band data. We do not find such an effect in the 
current data. Instead, the different spectrophotometric 
channels all have a similar level of noise relative to what 
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Fig. 4. — Raw spectrophotomctric light curves of WASP-19 (circles) and the corresponding correction functions (red lines) for the 
secondary eclipse observations obtained on UT 2012 April 5/6. The grey shaded region indicates the time where the planet is occulted 
(first to fourth contact). 
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is expected. 

3. ANALYSIS 

We fitted the spectrophotometric light curves for 
WASP-19 with transit and eclipse models multiplied by 
normalization and decorrelation (the quadratic function 
of time described in § 2.3) functions. We used th e exact 
analytic formulae given bv lMandel fc~Agol (pOOl for the 
transit and eclipse models. The parameterizations for 
these models are described in the subsections below. We 
assumed the planet is in a circular orbit, which is sup- 
ported by the existing radial velocity data and secondary 
eclipse times (lAnderson et al.ll2010t IGibson et all 120101; 
Hellier et al] 120111: IBurton et al.l 120121: lAnderson et al.l 
20131) . The orbital period of the planet (P) was fixed 
to the value determined from the transit times as de- 
scribed below. The normalization and decorrelation 
functions added an additional three parameters for each 
light curve. 

We identified the best-fit models and corresponding 
parameters for the analyses descri bed below using a 
non- linear least squares algorithm (Markwardt 2009). 
We determined the confidence intervals on the best- 
fit parameters using a residual permutation bootstrap 
( "prayer bead" ) algorithm. The uncertainties in the light 
curve points estimated during the reduction process were 
scaled to give a reduced % 2 = 1 for the best-fit model. The 
scaling factors range from 2.0 to 2.9. The uncertainties 
derived from the residual permutation bootstrap range 
from 1.0 to 2.5 times larger than the errors derived from 
a Markov Chain Monte Carlo technique, which suggests 
that correlated noise is the limiting factor for the obser- 
vations. 

3.1. Transits 

3.1.1. Model parameterization 

We parameterized the model for the transits using the 
square of the planet-to-star radius ratio ((Rp/R*) 2 ), the 
orbital inclination of the planet (i), the impact parameter 
(6 = a I R* cos i), quadratic limb darkening coefficients 
(71 and 72), and the central transit times (T c i). We 
used the same values for i and b for all the light curves. 
We assumed that the transit times were the same for all 
the spectrophotometric channels for a given event. The 
(Rp/R*) 2 values were determined for each spectrophoto- 
metric channel, and were assumed to be same for a given 
wavelength for both transit data sets when the two data 
sets were analyzed together. 

3.1.2. Limb darkening 

Determined transit depths from light curve fitting are 
strongly correlated with the adopted limb darkening de- 
scription even at the near-infrared wavelengths consid- 
ered in the current study. To bring additional constraints 
on the light curve fits and potentially increase the sen- 
sitivity of our analysis to the transmission spectrum of 
the planet's atmosphere, we estimated quadratic limb 
darkening coefficients using m odel atmospheres com- 
puted with the PHOENIX code (jHauschildt et al.lll999ft 
for the stellar parameters (T eff=5500K, log g — 4.5, 
and [M/H] =0.0) estimated bv lHebb et all (l2010j) which 
are consistent with the results from IDovle et al.1 (I2013T) . 
However, the light curve fits were significantly worse 



when using these coefficients as compared to when us- 
ing a single linear coefficient determined solely from the 
data themselves. For example, the models predict 71 
values ranging from 0.3 to 0.4 when 72 is set to zero, 
while the best-fit values with no constraints are below 
0.2 except in the bluest channel (see Tabled]). The dis- 
agreement between the theoretical and empirical esti- 
mates was not ameliorated when using model stellar at- 
mospheres with the WASP-19's nominal T e ff increased 
or decreased by 150 K. This is not an unusual situ- 
ation when modelin g high-precision light curves (e.g., 
iKnutson et al.l l2007t ). and likely stems from the inade- 
quacy of ID hydrost atic models to acc urately represent 
stellar atmospheres ()Havek et al.l 120121 ) . For all analy- 
ses presented in this paper we allow 71 to vary with no 
outside constraints, and we fix 72 = because our tests 
indicated that the data do not warrant an additional pa- 
rameter. 

3.1.3. Determining the transit parameters 

The first step in our analysis of the transit light curves 
was to estimate improved parameters for the system be- 
cause the MMIRS data are significantly more precise 
than existing data. We began by fitting the two tran- 
sit events separately with all the parameters free. The 
determined parameters were consistent between the two 
data sets. We then performed a combined analysis of 
the data sets requiring the parameters to be consistent 
for the events as described above. The determined phys- 
ical values for the transit in this case were also consistent 
with the most recently reported values in the literature 
(|Anderson et alj 120131 ) . To take advantage of the con- 
straints offered by the previously obtained data, we in- 
cluded priors for our final fit on i (79.42° ± 0.39°) and 
b (0.656 ± 0.015 ) based on the determined parameters in 
lAnderson et~aTl (|2013h . 

The individual transit times determined by fitting the 
MMIRS data are given in Tabled The determined tran- 
sit parameters, and a new ephemeris based on fitting the 
new transit times alon g with the previously pu b lished 
times from iHebb et al.1 (|201Ct ) and I Hellier et afl ()2011[ ) 
are given in TableO Note that we do not use the ad- 
ditional transit t ime for WASP- 19b av ailable in the lit- 
erature given by iDragomir et alj (|201lD for determining 
the revised ephemeris because it has a relatively low pre- 
cision (51s). None of the four considered transit times 
deviates from the new ephemeris by more than la. The 
average determined (Rp/R*) 2 = 0.0207 for the MMIRS 
data. 

3.1.4. Determining the transmission spectrum 

When assessing a planet's transmission spectrum from 
transit light curves, the relative depths are typically the 
key parameters to determine rather than the absolute 
depths because the pressure level probed by the data 
is often unknown, and thus comparison of the data to 
theoretical models for the planet's atmosphere must nec- 
essarily include a free offset in the overall level. We re-fit 
the MMIRS transit light curves after refining the transit 
parameters with i and b fixed to precisely estimate con- 
fidence intervals on the relative transit depths. We also 
fixed the transit times to the prediction of the ephemeris 
given in Table[2] because there is no evidence of transit 
timing variations in this system. The determined transit 



Bean et al. 



1.00 - 



0.98 - 



0.96 - 



J? 0.94 - 
> 

*| 0.92 1 

Q) 

0.90 h 



0.88 - 



0.86 




t — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — r 

0^ ^^f^jSfJI^^^^f^ - - 



V) 

o 

JO 
'(f) 
1) 
OH 



-4 
4(7" 

- 

-4 
4 

0h 
-4 

4 

0h 
-4 

4 

oh 

-4 
4 

0h 
-4 

4 

Oh 
-4 

417" 

- 

-4 
4 



_1_ 



_1_ 



_1_ 



_1_ 



- 1 — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — |- 

_l i i i i I i i i i I i i i i I i i i i I i i i i L 



- 1 — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — |- 

_l i i i i I i i i i I i i i i I i i i i I i i i i L 



-| — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — |- 
J i i i i I i i i i I i i i i I i i i i I i i i i L 



- 1 — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — |- 

_l i i i i I i i i i I i i i i I i i i i I i i i i L 



i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — c - i — i — i — |- 
i i i i I i i i i I i i i i I i i i i I i i i i L 



- 1 — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — i — i — i — |- 
J i i i i I i i i i I i i i i I i i i i I i i i i_L 



- 1 — i — i — | — i — i — i — i — r 

***** 

_i i i I i i i i L 



1 — i — i — i — i — | — i — i — i — i — p 



4 U I I I L±J I I ■ - f I I I I I I I I I I I I I I I I I I I L 



-3 



-2 



-1 



-3 -2 



-1 



Time from mid-transit (hours) 



Fig. 5. — Left Transit light curves (circles) and best-fit models (lines) for the WASP-f9b MMIRS data. The data for the two transits 
have been combined and binned to a sampling of one minute. Right Residuals from the best-fit models. 

with common F p /F* values. The F p /F* values deter- 
mined from this analysis are given in Tabled and the 
determined eclipse times are given in Table[3] The result- 
ing normalized and decorrclatcd secondary eclipse light 
curves with best-fit models and residuals are shown in 
FigureEl 



TABLE 4 

Transit Depths and Limb Darkening Coefficients 



Wavelength (/im) 



2 a 



f .25 - f .33 
1.40 - 1.50 
1.50 - 1.60 



1.60 
1.70 
1.95 
2.05 
2.15 
2.25 



1.70 
1.80 
2.05 
2.15 
2.25 
2.35 



0.02060 
0.02117 
0.02106 
0.02056 
0.02030 
0.02060 
0.02089 
0.02070 
0.02048 



2.8e-04 
3.1e-04 
1.5e-04 
3.0e-04 
1.6e-04 
3.1e-04 
1.0e-04 
1.8e-04 
4.7e-04 



0.28 
0.15 
0.15 
0.14 
0.12 
0.11 
0.16 
0.09 
0.07 



0.06 
0.07 
0.05 
0.05 
0.05 
0.07 
0.03 
0.04 
0.09 



a From an analysis with z, b, and 72 fixed to 78.73°, 0.681, 
and 0.0, respectively 

depths and limb darkening coefficients from this anal- 
ysis are given in Table|4] The resulting normalized and 
decorrelated transit light curves with best-fit models and 
residuals are shown in Figure[5J 

3.2. Secondary eclipses 

The free parameters for the secondary eclipse models 
were the planet-to-star flux ratios (Fp/F*) and central 
eclipse times (T C 2). We assumed the planet is a uniform 
disk. We used the values for i and b (see Table[2]) and 
the average value for the (R p /R*) 2 values (0.0207) deter- 
mined from the transit modeling. The eclipse times were 
the same for all the spectrophotometric channels for a 
given event. The F p /F* values were determined for each 
spectrophotometric channel. 

We first fit the light curves for the two secondary 
eclipses separately to check for systematic errors. The de- 
termined Fp/F* values for the spectrophotometric chan- 
nels were consistent between the different events. We 
then performed a combined analysis of the two data sets 



3.3. Possible influence of stellar activity 

Th e star WASP- 19 is known to be active (|Hebb et al.l 
I2010D . a nd this could potentia l ly infl uence our measure- 
ments. iTregloan-Reed et al.l ((2013) have recently re- 
ported the detection of star spot crossings in transits 
of WASP-19b that were observed in 2010. We do not 
see evidence for spot crossing in our transit light curves. 
However, the presence of spots on the host star could 
result in different transit and eclipse depths measured 
at different epochs even when they are unocculted be- 
cause the average brightness of the stellar disk might be 
changing. 

We performed some calculations to determine if our 
measurements could be significan tly influenced by unoc- 
culted spots. iHebb et alT (2010) reported photometric 
variability of up to 0.8% over the presumable 10.5 d ro- 
tation period of the star from the WASP-South discov- 
ery dat a (effective central w avelength of approximately 
550 nm, iPollacco et aLll2006fl . To estimate the maximum 
likely effect, we consider the case that the s pots leading 
to the photometric variability observed by IHebb et al.l 
(2010) are on the Earth- facing side of the star for one of 
the observations, and either not present or on the oppo- 
site side of t he star at the other epoch. Appl ying the for- 
malis m from lDesert et all ((201 If ) (Eq. 8 from iBerta et al.l 
l201ll) . and assuming the spot and star radiate like black- 
bodies with a temperature difference of 300 K and that 
the surface area fraction covered by spots is 4%, we find 
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TABLE 5 
Secondary Eclipse Depths 



Wavelength (fim) 



F P /F* 



Fig. 6. — Left Secondary eclipse light curves (circles) and best-fit models (lines) for the WASP-19b MMIRS data. The data for the two 
eclipses have been combined and binned to a sampling of one minute. Right Residuals from the best-fit models. 

from the mean transmission spectrum of much more than 
a few scale heights at this resolution. The standard devi- 
ation is 1.3 H, and the maximum spread between points is 
3.9 H. This suggests that low-resolution near-infrared ob- 
servations of planets like WASP-19b must achieve spec- 
trophotometric precisions on order of one atmospheric 
scale height or better to confidently detect molecular fea- 
tures. 

The second thing that can be seen from the data alone 
is that a flat line is a poor fit to the measurements. The 
best-fit flat line gives x 2 — 16-8 for 8 degrees of freedom, 
which has a 3% probability to happen by chance. The 
.ftf-band region is featureless within the limits of the pre- 
cision of the data. This is to be expected given that CH4 
is the main opacity source at these wavelengths and this 
molecule is not expected to be present in the atmosphere 
of such a hot planet like WASP- 19b assuming chemical 
equilibrium holds. However, the ii-band region of the 
spectrum exhibits variation. 

Theoretical models for WASP-19b's transmis- 
sion spectrum calcula t ed u sing the methods of 
Madhusu dhan &: Seagerl (|2009h are shown compared 
to the data in Figure[5] We considered two different 
solar metallicity models for comparison with the data: 
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0.00312 ± 9.1E-04 



that this situation would result in a 1 x 10 -4 change in 
the transit depth at 1.5 /im. The same situation would 
result in a change of 1 x 10~ 5 in the eclipse depths. These 
values are less than our measurement precisions. There- 
fore, we conclude that stellar variability likely has not 
influenced our measurements, and that we can combine 
our secondary eclipse data with previous data to put joint 
constraints on the properties of the planet's atmosphere. 

4. CONSTRAINTS ON THE PLANET'S ATMOSPHERE 

4.1. Transmission Spectrum 

The transit depths with relative errors (i.e. assuming 
fixed i and b, see §3.1.4 ) determined from the analyses of 
each transit separately and in combination are shown in 
Figure[7l The depths from the two transit observations 
all agree within 2.6a, and six of the nine values agree 
within l.ler. 

The transit depths for the combined analysis relative 
to the average value and in terms of the estimated scale 
height (H ) of the atmosphere (546 km) are shown in Fig- 
ure[H Two things can be seen from the data without 
comparison to models. First, the data rule out variations 



one with a roughly solar carbon-to-oxygen ratio of 0.5 
(|Asplund et al.|[2009t) . and one with a carbon-to-oxygen 
abundance ratio of 1.0 that was motivated by the 
recent idea that some h ot-Jupiters might ha v e such 
an abundance pa ttern (Madhusudha n et al.l I2011bt 
lMadhusudh1m1l2rjl2r ). 

The theoretical models are consistent with the data in 
the sense that variations of only a few scale heights are 
expected at this resolution. The solar composition model 
is actually a worse match to the data than the flat line. 
After adjusting the model by subtracting a fitted offset 



10 



Bean et al. 




1.2 1.4 1.6 1.8 2.0 
Wavelength (/iim) 



2.2 



2.4 



Fig. 7. — Transit depths with relative errors (i.e. assuming fixed 
i, b, and 72, see §3.1.4 ) determined from the analyses of each 
transit separately and in combination. The transit depths for the 
individual transits are shown offset in wavelength for clarity. 

to account for the unknown pressure level probed by the 
observations, the solar abundance model gives \ 2 = 24.3 
for 8 degrees of freedom. This poor fit suggests that the 
data are inconsistent with the solar abundance model at 
3.1(7 confidence. This inference hinges critically on the 
points at 1.45 and 1.75 /im. The main features in the 
solar abundance model are due to H2O, and these two 
channels should both show deeper transits compared to 
the two channels between them if this molecule is present 
in the planet's atmosphere. 

The carbon-rich model, while matching the data bet- 
ter than the solar abundance model, is only a marginally 
better fit than the flat line with % 2 = 15.4 after deter- 
mining a best-fit offset. The main spectral features in 
the carbon- rich model are due to HCN and H2O, though 
the H2O abundance is lower in the C-rich model than in 
the O-rich model by over a factor of 10. The molecule 
HCN has not been detected in an exoplanet atmosphere 
before, though several theoretical st udies have predicted 
its exi s tence in C-rich atmospheres ([M adhusudha n et al.1 
I2011at iKopparapu etalll20ll iMoses et al 1120131 ). Our 
observations favor the presence of this molecule, but the 
carbon-rich model does not provide a good global fit. 
Furthermore, even if a model could be found that pro- 
vided a good fit to the data, there is only a 2.1a confi- 
dence on the detection of spectral features. 

4.2. Emission Spectrum 

The secondary eclipse depths determined from the 
analyses of each eclipse separately and in combination 
are shown in Figure[5] The significance of the detections 
in the combined analysis range from 2.2a (the J- band 
point) to 14.4cr (the point at 2.1 /im). The data show 
the characteristic increase in the planet-to-star flux ratio 
with wavelength that is expected for hot- Jupiter spectra 
dBurrows et all 120051: iFortnev et all 120051: iSeager et all 
|2005[) . The depths determined from analyzing the sec- 
ondary eclipse events separately don't agree as well as 
the depths determined from analyzing the two primary 
transit events separately. This is likely due to the dif- 
ficulty of accurately estimating confidence intervals on 
low signal-to- noise detections. The eclipse depths from 
the two different events do all agree within 2.9a, and the 
combined analysis is likely more robust than the analyses 
of the individual events. 



The secondary eclipse depths from the combined anal- 
ysis of the two observed events are shown in FigurelTOl 
along with previously published data and example theo- 
retical models. Previous secondary eclipse measurements 
of WASP-19b using photometric techniques have bee n 
made with with IRA C on Syitzer (lAnderson et al.112 013). 
HAW KI on the VLT (lAnderson et a lJl2010 UGibson et al.l 
20101: iLendfeT^ I20T1, ULTRACAM on the NTT 



Burton et al. (12013 ). EulerCam on the Euler-Swiss tele- 



scope (|Lendl et al.ll2013D . and TRAPPIST ()Lendl et al.l 
I2013D . Our data overlap broad-band H a nd narrow-band 
K mea sure ments previous l y rep orted by lAnderson et al.l 
(|2010D and iGibson et al.l (|2010l ). respectively. We per- 
formed an analysis of our data with the bandpasses set 
to closely match the wavelengths of those previous mea- 
surements to check for consistency. We find that the 
results are consistent within 1.7 and 1.3cr for the H- and 
K-band points, respectively. 

The combined data set presented in FigurelTUl provides 
coverage of the planet's spectral energy distribution from 
0.9 to 8.0 /xm that is matched only for a few other exo- 
planets. We performed a spectral retrieval analysis on 
the combinatio n of our data with the IRAC po ints using 
the methods of IMadhusudhan fc Seagerl (|2009fl to inves- 
tigate what constraints can be placed on the composition 
and structure of the planet's atmosphere. 

An acceptable fit to the data can be obtained assuming 
an isothermal atmosphere with a temperature of about 
2250K (x 2 = 19.2 for 12 degrees of freedom). While a 
perfectly isothermal atmosphere over the entire day-side 
of the planet may be unrealistic, an effectively isothermal 
1-D temperature-pressure profile could arise if the planet 
has a thermal inversion where heating due to absorption 
of incident radiation from the host star causes the tem- 
perature to increase at lower pressures ins t ead of decrease 
(IHubenv et al.H2003t iBurrows et al1l2007t IKnutson et al.l 
I2008D . The planet's atmosphere could appear isother- 
mal because the heating yields only a modest change in 
temperature over the pressures probed by the observa- 
tions. In this case, there could be a large temperature 
increase at lower pressures that aren't being probed given 
the resolution and precision of our data. An effectively 
isothermal atmosphere with a temperature of roughly 
3000 K has also recently been sugg ested for the planet 
WASP-12b (jCrossfield et al.l l2012af ). Higher resolution 
and higher precision data are needed in both cases to de- 
tect spectral features and further constrain the planets' 
temperature-pressure profiles. 

The data for WASP-19b disfavor a thermal inversion 
where there is a large temperature increase to lower pres- 
sures in the observed part of the atmosphere in both 
oxygen-rich (C/O = 0.5, see Figure [TU]) and carbon-rich 
(C/O = 1.0) c omposition mod e ls. Th is is consistent with 
the finding of lAnderson et a l. (2013). In the presence of 
a strong inversion in the observable atmosphere, absorp- 
tion lines due to H2O and CO that are prominent in the 
IRAC bands would reverse to emission and yield a poor 
fit to the Spitzer data. 

We also explored models for the planet's atmosphere 
with no thermal inversions, i.e. where the tempera- 
ture decreases monotonically with pressure. Figure lTUl 
shows two example models with different chemical com- 
positions: one with an oxygen-rich solar composition 
(C/O = 0.5) and another with a carbon-rich composi- 



11 




Fig. 8. — The transmission spectrum of WASP-19b in terms of the relative transit depth (left axis) and scale height of the planet's 
atmosphere (right axis). The measureme nts are given as the black circl es. Also shown are models (lines) for the planet's transmission 
spectrum calculated using the methods of Madhusudhan & Scagcr (2009). One model has a solar carbon-to-oxygen ratio ("O-rich"; solid 
red line) and the other model has a carbon-to-oxygen ratio of 1.0 ("C-rich"; dashed blue line). The colored points give the model values 
binned over the bandpass of the observations. The models have been adjusted up and down to give the best-fit to the data. 
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Fig. 9. — Secondary eclipse depths determined from the analyses 
of each eclipse separately and in combination. The depths for the 
individual eclipses are shown offset in wavelength for clarity. 

tion (C/O = 1.0). The C-rich model provides a better 
fit to the data (x 2 — 15.9) compared to the O-rich model 
(x 2 = 30.7), although the significance is low due to the 
low number of degrees of freedom, which is only three 
(the models have 10 free parameters). The absorption 
features of H 2 that are expected in the O-rich scenario 
are not obvious in the MMIRS data, which is consis- 
tent with what is seen in the transmission spectrum. On 
the other hand, though the C-rich model provides an 
acceptable fit to the data, it is not much better than an 
isothermal model with any composition given the current 
precision of the data. 

The isothermal and carbon-rich non-inversion models 
fitted to the IRAC and MMIRS data also provide a rea- 
sonable fit to the data at other wavelengths that were not 



includ ed i n the fit. The z'-ba nd points of iBurton et al.1 
( 2012) and lLendl et al.l (|2013f ) are inconsistent with each 
other at 2 Aa. The fit t ed m odels happen to go right 
through the lLendl et al.l (120131) p oint. The deeper eclipse 
observed by IBurton et al.1 (|2012f ) suggests a higher tem- 
perature that is difficult to reconcile with our data. 

5. DISCUSSION 

We have presented the first near-infrared transmis- 
sion and emission spectroscopy measurements of a hot- 
Jupiter obtained using the technique of multi-object 
spectroscopy with wide slits. The data do not yield 
strong detections of spectral features, yet they do provide 
some constraints on the physical properties of WASP- 
19b's atmosphere. We rule out broad spectral features 
in the planet's near-infrared transmission spectrum that 
would arise from probing more than a few scale heights in 
atmospheric pressure. The combination of our thermal 
emission measurements with extant Spitzer IRAC data 
suggests that the planet does not have a thermal inver- 
sion yielding a large increase in temperature to lower 
pressures in the observed part of its atmosphere. The 
emission data are consistent with a model for the planet's 
atmosphere that is isothermal over the observed part of 
the atmosphere, and with a model that has decreasing 
temperature with pressure and a C/O value that is sub- 
stantially larger than the solar value. 

A key question raised by our observations is why 
isn't the obtained precision higher than it is? The cur- 
rent data, and our previous near-infrared observations of 
GJ 1214b (jBean et al.ll201lD . fall short of delivering light 
curve residuals better than a factor of two to three times 
the photon-limited expectations. Furthermore, corre- 
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Fig. 10. — The thermal emission spectrum of WASP- 19b in terms of the planet- to-star flux ratio (seconda ry eclipse depth ) . Th e 
measurement s pres e nted in this paper a r e shown as the black circle s. Previous repor ted measurements from [Anderson ct al. (2010), 
I Gibson et"aH 12010). Burt oneFaT] 1 120121 ). [Anderson et all 1 12013 ). and lLendl et al.l j2013) are also shown as the blue {Spitzer IRAC) and 
purple (ground-based) circles. Th e green and red lines repr esent models without thermal inversions that were fit to our data and the Spitzer 
IRAC data using the methods of Madhusudhan & Scaecr (2009). The green line is a roughly solar composition model (C/O = 0.5), and 
the red line is a carbon-rich model (C/O = 1.0). The grey line shows an example model with a thermal inversion and solar composition. 
The green, red, and grey points show the respective values of the models integrated over the bandpasses of the photometric points. The 
integrated model values for our spectroscopic data are not shown for clarity. The inset shows the corresponding temperature-pressure 
profiles for the models. The dashed line is a model of the planet assuming it radiates as a blackbody with a temperature of 2250 K. The 
dotted lines are blackbody models for the planet with T= 1800 (lower line) and 2900 (upper line) K. 



lated noise in the light curves degrades the precision of 
the transit and eclipse depth measurements even more. 
This is in contrast with observations that we have done 
at optical wavelengths where we have obtained residuals 
as small as a few hundred ppm per minute and that are 
within a f ew tens of percent of the photon-limited preci- 
sion (e.g., IBean et all 120101 ) . Our study of the MMIRS 
detector non-linearity described in § 2.2.2 yielded surpris- 
ing results, and this leads us to believe that the quality 
of near-infrared detectors could be responsible for the 
lower than expected precision being obtained at these 
wavelengths. However, we can't speculate as to what the 
underlying physical reason might be. We do note that 
if these observations represent the limit obtainable with 
ground-based near-infrared spectroscopy, then roughly 
six transits would have to be observed to detect spec- 
tral features in the transmission spectrum of a planet 
like WASP- 19b at better than 5cr confidence. 

The near-infrared is clearly an important wavelength 
region for exoplanet atmosphere measurements because 



of the presence of many different molecular bands. How- 
ever, ground-based near-infrared spectroscopy to mea- 
sure these molecules is challenging due to the presence 
of water vapor in Earth's atmosphere. Numerous and 
strong telluric water lines sculpt the near-infrared spec- 
trum and limit observations to the canonical Y-, J-, H-, 
and if -bands. Observations at the wavelengths between 
these bands are simply not possible from the ground 
because no light from astronomical sources reaches the 
telescope. The presence of telluric water vapor further 
complicates matters because it produces lines that even 
contaminate the windows where light can get through. 
These lines are expected to be strongly variable over 
timescales relevant for transit observations due to chang- 
ing conditions in the atmosphere above an observatory. 
This is borne out by our observations. We clearly see 
stronger variations in the atmospheric transparency as 
measured by the reference stars between the edges and 
centers of the H- and AT-bands (see Figured]). 
Variable telluric water vapor has been identified as a 
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possible reason for the strong feature seen in the emis- 
sion sp ectrum of t he planet HP 18 9733b that lSwain et al.l 
( 2010) presented (jMandell et all 120111 ). The technique 
of observing only the transiting planet system, which 
ISwain et al.l (|2010f ) use, requires correctin g for telluric 
variat ions using the data themselves. The ISwain et al.1 
(2010) approach is to filter out these variations based on 
the assumption that they are correlated with wavelength 
and time. In contrast, the multi-object technique we use 
enables corrections for variations in Earth's atmospheric 
transparency at every wavelength for every exposure us- 
ing an external reference. Therefore, our data should 
be insensitive to the problem of variability in water va- 
por lines assuming the variations are common mode for 
sources within a few arcminutes of each other on the sky. 

The measurements we have presented here represent 
one of the few cases of ground-based exoplanet atmo- 
sphere observations with repeated observations. The 
moderate level of disagreement seen between the results 
from the separate analyses of the different events illus- 
trates the potential pitfalls of making statements about 
the nature of exoplanetary atmospheres at even a 3er for- 
mal confidence level when the inference critically depends 
on one or two points that are derived from observations of 



a single event. Our data require little systematic decor- 
relation and we used a well-regarded algorithm for esti- 
mating confidence intervals in the face of correlated noise 
(residual permutation bootstrap), but yet we still only 
see repeatability in our transit and eclipse depths at the 
2-3cr level. One strength of ground-based measurements 
is that the observations are easier to repeat due to the 
generally lower time pressure on the telescopes compared 
to space-based facilities. We want to encourage observers 
pursuing ground-based observations to take advantage of 
this opportunity, and we also want to encourage telescope 
time allocation committees to look favorably on proposed 
repeated observations not just because they would have 
a higher formal signal-to-noise, but also because the re- 
sults will likely be more robust. 
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ministered by the NASA Exoplanet Science Institute 
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at Las Campanas Observatory. 
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